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1. INTRODUCTION 

The direct conversion transceiver due to its simplified and low-cost structure is the most popular 
architecture and the appropriate choice often used in wireless communication systems [1]. Despite the 
simplicity and integrability of this structure, there are problems due to its practical implementation of analog 
components especially in the applications of advanced modulation schemes with high peak to average power ratio 
(PAPR). The nonlinear nature of the power amplifier on one hand and the in-phase and quadrature-phase 
asymmetry impairments of the modulator, on the other hand, cause in-band signal distortions and adjacent 
channel interference and dramatically reduce the performance of the transceiver system [2]. Moreover, the 
advanced signal waveforms such as orthogonal frequency division multiplexing (OFDM), with wide bandwidth 
and increased PAPR, cause the negative effects of power amplifier nonlinearity to be worse. Therefore, using 
linearization techniques to improve the efficiency of the transmission signal is necessary [3]-[5]. 

The digital pre-distortion (DPD) technique as one of the most effective methods, is typically used to 
overcome the non-idealities. To utilize this method, it is necessary to model the nonlinear behavior of the 
system with a proper mathematical model such as the memory polynomials [6]—[8], the Wiener Hammerstein 
model [9], Volterra series [10], [11], the optimized structures with modified algorithms [12], [13] and 
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orthogonal polynomial [14], [15], which consider the memory effects along with the static nonlinearity. 
Among the many models that have been proposed in the literature, the memory polynomial model is a true 
parametric model which represents a simplified case of the Volterra model as only diagonal terms are 
considered and the others are zero. To estimate the model coefficients, the least-squares (LS) approximation 
is one of the most widely used approaches. In this method, the inversion of an observation matrix must be 
calculated, which leads to arise an instability problem, dependent on the size of the matrix and the statistics 
of the input signal. To prevent this instability problem, Raich and Zhou proposed a set of orthogonal basis 
functions based on the Laguerre polynomials in [14]. The proposed model worked well only for the signals 
with Rayleigh distribution. But depending on the signal distribution, the conditioning of the observation 
matrix drastically deteriorates and decreases the DPD performance. As nonlinear order and memory depth 
increase, the problem becomes more serious. Aziz et al. in [16] proposed a rational-function—based model to 
jointly attenuate the power amplifier (PA) nonlinearity and in-phase/quadrature-phase (IQ) impairments. The 
proposed model improved the numerical instability problem of DPD by a reduction in the number of 
coefficients and so reduced the dispersion coefficients and condition number. In [17] a distributed model 
using conventional memory polynomials was presented. The proposed model used a two-block structure for 
DPD, the first block was a memory polynomial model to reduce the nonlinear effects of the PA, and the 
second one was a mild nonlinear model with memory to compensate the modulator's impairments. It was 
claimed that by distributing nonlinearity and memory depth, a considerable reduction in matrix conditioning 
and coefficients dispersion was obtained. The obtained numerical stability due to the reduction in the number 
of pre-distortion coefficients in [16], [17], is very fragile and as the nonlinear order of the system intensify, 
the number of pre-distortion coefficients inevitably increase, resulting in the changes in the matrix 
conditioning and numerical instability arise. An orthogonal polynomial model was introduced by Wang et al. 
in [18] by using the combined least-mean-square (LMS) and recursive least square (RLS) algorithm to model 
the nonlinearity of PA with memory effect. Although this method improved the DPD performance, it was 
limited to the input signal values distributed in [0, 1], and the combined LMS and RLS algorithm increased 
the computational complexity. Manai et al. [19] proposed a numerically stable digital pre-distorter based on 
the Gegenbauer polynomials to linearize the radio frequency (RF) power amplifiers without the IQ 
impairments compensation. Although the Gegenbauer polynomials were orthogonal and the pre-distorter 
based on the Gegenbauer basis functions was numerically stable, it worked only on the real-valued variable 
on the input values of the interval [—1, 1] and only considered the PA nonlinearity while all the other 
impairments of the transmitter affected on the DPD performance. 

In this paper, a set of new orthogonal basis functions is introduced based on the Chebyshev 
expansion to joint compensate for the non-linearity and the memory effects of the power amplifier as well as 
frequency-dependent impairments of the direct conversion transmitters. The new DPD structure is intended 
to solve the numerical instability problem that may arise through the conventional polynomial model 
approximation, without any limitation on the input signal statistics. Therefore, by increasing the nonlinear 
order of the system and as a result increasing the number of pre-distortion coefficients, the value of the 
conditional number remains small and numerical stability is guaranteed. We achieved the best numerical 
stability merits with the least numbers of coefficients and computational complexity in comparison to the 
references’ literature. 

This paper is presented in five sections: section 2, describes the mathematical model of PA 
nonlinearity, then proposes the Chebyshev pre-distortion method to alleviate all the transmitter imperfections. 
In section 3, the numerical stability and performance evaluation of this work are analyzed. The simulation 
results demonstrate the operation of this work in section 4 and finally, the conclusions are given in section 5. 


2. MODEL DESCRIPTION 

In this section, the memory polynomial (MP) model is presented to model the nonlinearity of PA 
with memory effect, then the numerical instability problem of the MP model for DPD is explained. The new 
DPD model based on Chebyshev polynomials is proposed and leads to making the new joint pre-distorter to 
compensate for the PA nonlinearity and IQ impairments. 


2.1. Modeling nonlinearity of power amplifiers regarding the memory effects 

Power amplifier as a main element of the transmitter chain has a nonlinear characteristic and distorts 
the transmission signal which results in spectral regrowth in digitally modulated signals. This phenomenon 
causes signal interference in neighboring channels. To correct the destructive effects of power amplifier 
nonlinearity, the behavior of this component must be mathematically modeled. To have a more complete 
model, the effects of memory must also be considered. For this reason, a memory polynomial model is 
derived by (1): 
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y(n) = Dono Dono C2p+1,4 P(n — q) (n — q)? (1) 


where v(n) is the complex base-band signal at the input of PA, P and Q represent the polynomial order and 
memory depth of the model respectively, and y(n) is nonlinear PA output. The complex-valued coefficients 
C, are given as follow, which are extracted from an actual Class AB PA [20]. 


C19 = 1.0513 + 0.0904; C11 = —0.0680 — 0.0023; cz = 0.0289 — 0.0054; 
C39 = —0.0542—0.2900j c34 = 0.2234 + 0.2317 C39 = —0.0621 + 0.0932; 
cso = —0.9657 — 0.70287 Cay = —0.2451 — 0.3735f Cap = 0.1229 + 0.1508 


2.2. Numerical instability problem in the model identification 
The conventional polynomial model of a non-linear system can be expressed as in (2): 


y(t) = DK, belel x(t) (2) 


where x(t) and y(t) are the baseband input and output of the nonlinear system and b, is the unknown complex 


model coefficient. By replacing p(x) = |x|*~1x, the (2) can be rewritten as (3): 


y(t) = Ek=1 debe (x(t)) (3) 


where (x) is the conventional polynomial basis function. To extract the model coefficients b = 
[b,, b2, ..., bg]”, N samples of the measured input and output envelopes can be used as N-by-1 input/output 
data vector:x = [x(t,),...,x(ty)]", Y = [y(t,),..., y(ty)]", then define the following vector notation. 


P(X) = [be (x(t1)),-- PAND D= G2(x) + Oe] (4) 
The relation (3) can be written as Y = ®b and the least square estimate of the coefficient vector b is 
brs = (b' p tory (5) 


where (b' @ )~1@# is the Moore-Penrose pseudo-inverse of ® and ®# is its Hermitian transpose [14]. The 
inversion of the K-by-K matrix ®"@ can experience a numerical instability problem because the matrix 
©" is often ill-conditioned [14]. To overcome this problem, various methods have been proposed to 
convert the conventional polynomials to a set of orthogonal basis functions and thus improve the condition 
number of the observation matrix. But these orthogonal polynomials have some limitations in signal 
statistics. On the other hand, some classic orthogonal polynomials have very useful properties in solving 
mathematical and physical problems but are less commonly used in DPD structures. In the next section, we 
propose a set of orthogonal basis functions based on Chebyshev expansion to attenuate the problem of 
numerical instability arising during the conventional model identification method. 


2.3. Chebyshev polynomials as the pre-distorter model 

Digital pre-distortion based on memory polynomial model is one of the popular and efficient models 
that is used to compensate the PA nonlinearity but can experience the numerical instability problem that 
reduces the performance of the DPD. To overcome this problem, various methods have been proposed to 
convert the conventional polynomials to a set of orthogonal basis functions [14], [16]-{19] and thus improve 
the condition number of the observation matrix. But these orthogonal polynomials have some limitations in 
signal statistics or some of them are complicated in the calculation [14]. The orthogonal polynomials are a 
category of polynomials that satisfy the orthogonality relation (the inner product between two terms of the 
sequence must be zero.) on a specified interval [18]. In this section, an orthogonal polynomial model based 
on Chebyshev expansion is proposed for the digital pre-distorter modeling. The Chebyshev polynomials of 
the first kind are defined as in (6): 


To(x) = 1, T,(x) = x, 
Tns) = 2xT, (x) — Tn- (x) (6) 


where Ta(x), n = 0,1,2, ..., is the n™-order term of the polynomial. The original first kind Chebyshev 
polynomials formula is given as in (7) [21]: 
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T(x) = p P21) Loe pp-2m-1yp-2m, p = 1,2,3, 0 (7) 


m!(p—2m)! 


where |p/2] denotes the largest integer number less than or equal to p/2 and p is the order of nonlinearity. 


As stated earlier, the Chebyshev basis functions are orthogonal and their orthogonality is limited to 
the real-valued variable in [-1, +1]. Whereas baseband signals used in communication systems are typically 
complex with an interval [a, b] where b > a. So, they cannot be directly used as the basis functions in digital 
pre-distortion structure. Furthermore, the zero-order term in the Chebyshev polynomial model does not 
conform to the baseband model of a nonlinear system [14]. To make the proposed Chebyshev polynomial 
appropriate to any finite range of the communication signals (ie.a<z<b) and apply this orthogonal 


polynomial to the complex signals as a pre-distorter model, we propose the new variable Z by applying the 
linear transformation as 7 £ 2(r — a/b — a) — 1 with r = |z| and a = min(r),b = max(r), Then replace 
thex?~*™ in (7) by (2)?-*™z, to define T,,,(z) as in (8): 


i) =P ey er Oe for po dza (8) 


m!(p—2m)! 


where T,(z) = z; z is a complex variable of the input signal and is defined as z = x + iy; xy E€ R. 

Now, equation (8) defines the modified Chebyshev polynomial which can be used as a nonlinear 
model to design a new robust pre-distorter in the next section. The closed-form relation of the modified 
complex Chebyshev polynomial with the memory effects can be written as in (9): 


Ta @@) = Pigei do. CUM ae? Ce =a)? = 4), 
for p=1,2,3,... 
Ty g(2(t)) = Deg Z(t — q) (9) 


where Q is the memory depth and p is the highest non-linearity order. 


2.4. Proposed digital pre-distorter model of PA and IQ modulator 

Figure 1 shows the baseband model of a direct conversion transmitter with a joint digital 
pre-distorter. The schematic of the joint pre-distorter is shown with more details in the mathematic relations. 
It is based on the Chebyshev polynomial model which can linearize several different nonlinear models with 
memory. Based on the pre-distorter model in Figure 1, the mismatch compensator block consists of two finite 
impulse response (FIR) filters that compensate for the frequency-dependent IQ impairments of the modulator 
[22]. For the joint pre-distorter, we must combine the PA pre-distorter with the mismatch compensator and 
build a unit structure from both of them. The pre-distorter output is represented as in (10): 


up(n) = E£ o a7 (n).T,(s(n)) (10) 


f ' ' ; T T 
where Ty(s(n)) = [T14 (S0), Taa (S0), -s Topsaq(S@))] > aq = [ara n), a3q(N), a2p+14 (0| 
and q indicates the memory size. 

By converting the pre-distorter output in (10) to a vector form, we can show that equation as a 
closed-form in the following: 


u, =a" (n).T'(s(n)) (11) 


To take into account the modulator mismatches, an obvious form for pre-distorter of the modulator 
can be obtained by the following signal transformation [22]: 


x(n) = Gy(n)  up(n) + G2(n) @ up (n) +d (12) 
By using (11) in (13), x(n) is rewritten as in (13): 


x(n) = G (up (n)) + G2(up* (n)) +d l l 
= G;(up(n)) + Go(da. Up(N)) + G;(up(n)) + Go (df.uz(n)) +d (13) 
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where da and df are constant-coefficient related to the gain and phase mismatch of the modulator. So, after 
simplification we have: 


x(n) = G;(up(n) + us (n)) + Go(da. up(n) + dB.us(n)) + d (14) 
Gi(n) = Dio tex(n — k), Gon) = Lio bi x(n — D). (15) 
Suppose to L = K then x(n) can be rewritten as: 
K 
x(n) = Yd.) a" (n-HT'(s(n-W) +a" -0T s-k) | +d (16) 
k=0 


where d;,(n) denotes the FIR filter coefficients vector and its order are equal to the number of FIR filter taps. 
So, the final joint pre-distortion model formulation is written as (17). 


x(n) = 2o alq (n). Tka (S(N)) + pies ali, (n).T “zq(s(n)) +d (17) 


demodulator 


Figure 1. Baseband equivalent model of direct conversion transmitter employing joint pre-distorter 


Figure 2 has demonstrated the proposed model. It should be noted that the frequency-dependent //Q 
impairments are modeled as a linear time-invariant (LTI) system and the memory polynomial PA can be 
displayed as a two-box structure that contains an LTI system followed by a memory less nonlinear system. 
To build a joint model, we unified the LTI system of the I/Q modulator with the LTI portion of the power 
amplifier and obtained one merged LTI system which is connected to a memory less nonlinear system in the 
joint compensator. After this simplification, we have the final coefficients vector in order of 2(K + Q + 
1)(P + 1) +1, For example, if K = 1, Q = 2 and P = 3, then the order of joint pre-distorter coefficients 
will be 33. 


O+K 
D Aig (2)-Ty, (s()) 
kq=0 


Figure 2. The final schematic of the simplified joint digital pre-distortion model 


Joint digital pre-distortion model based on chebyshev expansion (Elham Majdinasab) 


3786 O ISSN: 2088-8708 


Joint DPD model coefficients are adaptively corrected using the indirect learning architecture (ILA) 
[23]. According to the architecture in Figure 1, the coefficients that are estimated for the pre-distorter and 
compensator training blocks in the feedback direction are equal to those of the pre-distorter and mismatch 
compensator blocks in the forward direction, therefore the discrete-time output signal of the training block 
can be written as: 


z(n) = bj, (MT; (w(n)) (18) 


where, b, (n) =[b,,"(7).b,,"(a),d! TL Gun) =T (wn), T) (w(n)) AY . 
By using the least square cost function, the training block output signal z(n) and joint DPD output signal 
x(n) will be close together to minimize the differential error: 
n 2 n 2 
J(n)= Ja D- = a bw = bj, (Tw) (19) 
l=1 


l=1 


where A is named the forgetting factor and its value is a real number in the domain 0</<1 (in our calculations 
A = 0.988). In the following, an RLS algorithm has been described [23] to minimize J(n) and estimate the 
optimal b,g(n). After finding the by,(n) at each iteration, they are replaced by akq (n) in (13) as the 
pre-distorter coefficients, and the process will be continued until x(n) be equal to the z(n), see Figure 1. 


2.5. RLS algorithm for parameters estimation 
The RLS algorithm is defined as following steps: 
— Initialization steps: 


QO) =6'T , b (0) =[10,...,0] 


Where J is a small positive constant for initialization (in our work ô = 0.008), I denote the identity matrix. 
— N” iteration: 
E A'Q(n-1) Ti," (W) 
1+ 47T (WDN -DT wn) 
e(n) = x(n) -by (n- DT; (w(n)) 
by, (n) =by(n-1)+k(n)e" (n) 
Qin) =A 'Q(n-1) -Ak T; (WAN -I) 


k(n) 


3. STABILITY ANALYSIS AND MODEL PERFORMANCE 

In this section, at first, the numerical stability and complexity of this method are analyzed. The 
condition number and the dispersion coefficient are two numerical analysis metrics to evaluate the numerical 
stability of the proposed model. Then the performance evaluation of the proposed scheme is demonstrated to 
validate the presented model. 


3.1. Numerical stability and complexity analysis 

As shown in section 2.2, the coefficients vector of the joint DPD model (bkq) can be calculated 
using the Moore-Penrose pseudo-inverse. This pseudo-inverse calculation is significantly influenced by 
changes in matrix conditioning [24]. The definition of the condition number depends on the norm type, since 
we used the RLS algorithm based on the norm-2 error cost function, the norm-2 condition number, K(A), is 
utilized to measure the numerical stability of the proposed model (A is the matrix to be inverted). Therefore, 
it is defined as the ratio between the highest (01,4) and the smallest (02,4) eigenvalues calculated for the 
Vandermonde matrix using single value decomposition [24]: 


(A) = ||AlloI|A“lp = 24 (20) 


02,A 
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In our calculations, T-*4 is the Vandermonde matrix includes the augmented power of the input 
signal. Figure 3 shows the condition number of the observation matrix vs. the polynomial order for the 
proposed Chebyshev DPD model in comparison MP model. As can be seen in the figure, by using the 
pre-distorter model, we have achieved a significant reduction in the matrix conditioning (A) for the LTE 
16 QAM signal. 

Another numerical analysis metric is the dispersion coefficient. The dispersion coefficient specifies 
the dispersion of the extracted coefficients vector over the domain. Higher dispersion coefficients represent a 
larger number of bits used to cover the entire range of coefficients in the digital signal processor (DSP) chips 
[24]. Figures 4 and 5 show the coefficient dispersion and the condition number of the observation matrix for 
each sample of the input signal during the ILA iterations. As shown in the plot, for the OFDM 16 QAM 
signal, the condition number has improved more than 80 dB in comparison to the MP model, and the 


dispersion coefficient parameter has a reduction of around 10-20 dB, which causes a considerable reduction 
in complexity. 
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Figure 3. Condition number of the different polynomials’ observation matrices for the LTE 16 QAM signal 


Coefficient Dispersion (dB) 


% 0.5 1 1.5 2 2.5 3 3.5 
Samples (N) x10 


Figure 4. Coefficient dispersion of the different polynomials’ observation matrices for the OFDM 16 QAM 
signal 


Condition Number(dB) 


0.5 1 1.5 2 2.5 3 3.5 
Samples(n) x10“ 


Figure 5. Condition number of the different polynomials’ observation matrices for the OFDM 16 QAM 
signal 
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The figures of merit for the complexity evaluation of the proposed joint DPD are the number of 
complex-valued parameters and the number of floating-point operations (FLOPs) required to implement the 
model. In this work, a new orthogonal Chebyshev DPD has been defined which its polynomial basis 
functions with odd-order nonlinearity are Ty4(s(1)) = [Ti ka (S0), T3,heq (S(O): Ts ra (S0), T; ka (S(M)), ale 
kq = 0,1,2,3. To calculate the number of required FLOPs, we have driven two relations based on the method 
used in [17, 25] for the Chebyshev DPD (based on the recursive equation in (2)) and the MP DPD (based on 
equation (1)). According to the results shown in Table 1, for a joint DPD with 33 complex-valued 
coefficients (K = 1, Q =2 and P = 3), we require 356 FLOPs at every step of the computation to 
implement the model. Although more FLOPs are required for Chebyshev DPD implementation in 
comparison to MP DPD, the proposed DPD builds strong numerical stability with the least coefficients and 
complexity. 


3.2. Model performance metrics 
To evaluate the linearity performance of the system, the normalized mean square error (NMSE) is 
commonly used [25], [26], which is formulated as follows: 


N a Z 2 
NMSE(dB) = 10 logio (r= aesirea Cu Y measured (n)| ) (21) 


Eral desirea (n)|? 


where the desired signal is denoted the input signal and the measured one is the PA output after the joint 
digital pre-distortion. The results are shown in Table 1, where the NMSE performance of the transmitted 
signal after PA linearization and compensation of IQ impairments for the proposed model has improved more 
than 3 dB in comparison to the MP model. 

A principal criterion to compute the power leakage due to the spectral regrowth into the adjacent 
channel is adjacent channel power ration (ACPR) which is given by [17]. 


(22) 


adj.channe Sf) 
ACPR ap) = —10logio (uem N 


Sinbana Sf) 


In the above statement, S(f) defines the power spectral density of signal, adj. channel, and in-band defines the 
adjacent channel and desired channel, respectively [17]. PA nonlinearity disturbance and I/Q impairments 
cause excessive distortion of the transmitter output signal, resulting in poor ACPR. However, by performing 
the proposed pre-distorter, the transmitter impairments are compensated and significant improvement in 
ACPR, about 16 dB, is achieved. For the input symbol Sin and the transmitted symbols Sry, the error vector 
magnitude (EVM) between the two is defined as in (23). 


Wee lSin®-Srx OI? 


1 F 
SlalSra(i)l? 


EVM = 


x 100% (23) 


The results are presented in Table 1 and after using the proposed DPD method for the class-AB PA, great 
improvement for EVM has been obtained. 


Table 1. Performance evaluation of Chebyshev DPD in comparison to memory polynomial DPD 


DPD No. of No. of Flops ACPR EVM% EVM% NMSE 
Model Coefficients (dB) (before) (after) (dB) 
Chebyshev 33 (K+Q+1)(2Px4+2+8(P+1))+8(K+Q 43.2 45.49% 0.66% 3034 

+1)(P +1) —4 = 356 dB 
104K +Q+1)(P+1)+P(K+Q+1)+8(K+Q -33.27 
MP 33 41)(P + 1) — 4 = 296 -41 45.49 % 1.3 % JB 


4. RESULTS AND DISCUSSION 

To demonstrate the performance of the common digital pre-distorter model, a signal transmission 
system is simulated in MATLAB to model a direct conversion transmitter with impairments. The system is 
excited by an OFDM 16-QAM single carrier modulation signal with about 4 MHz bandwidth and 7.68 MHz 
sampling frequency in the digital front-end. A memory polynomial nonlinear model of PA from [20] is 
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simulated which is related to a real class AB power amplifier. An indirect learning architecture is used to 
model digital pre-distorter structure. The values of the gain and phase imbalance used in the simulations are 
g=1.07 and ọ = 5°, respectively, and a dc offset which is equal to c=0.05+0.01i. The frequency-dependent 
mismatch is modeled by two filters as hi=[0.99 -0.1]" and hg=[0.98 -0.07]', for I and Q branches, 
respectively. The Chebyshev pre-distorter is employed as three/two order (nonlinear order/memory depth), 
with filter lengths of 2. We used 30,720 samples in the RLS algorithm to estimate the digital pre-distorter 
coefficients. Figure 6(a) shows the power spectral density (PSD) of the output signal after applying the 
proposed pre-distorter. It can be seen that without any pre-distortion, spectral regrowth will occur and causes 
high distortion in co-channel and adjacent-channel which eventuate a poor ACPR. According to the results, it 
is clear that by using the proposed pre-distorter, the ACPR has enhanced greatly, and the PSD of the PA’s 
output fully complies with that one which is from the input signal. These simulations are executed for the 
OFDM 16 QAM input signal and an evaluation of the performance is presented. It can be observed that the 
proposed joint pre-distortion model was able to compensate for the memory effects and nonlinearity of PA as 
well as IQ imbalance of modulator in direct-conversion transmitter (DCT). Based on the results, by applying 
this DPD structure, for high PAPR input signal (about 10.3 dB), in presence of transmitter non-idealities, 
ACPR more than 43 dB can be obtained which is about 3 dB more than the MP DPD model. When IQ 
impairments effects on input signal are considered, the PA output has a more complicated nonlinearity, 
which can be seen as amplitude and phase modulation, known by the amplitude modulation to amplitude 
modulation (AM/AM) and amplitude modulation to phase modulation (AM/PM) parameters, respectively. 
The characteristics of the transmitter are shown in Figure 6(b). After linearization, the desired AM/AM and 
AM/PM characteristics have been obtained with the least model coefficients. In Table 2, we have compared 
the proposed method with the methods presented in the references to show more the superiority and 
efficiency of this method. All the table parameters are obtained from the values reported in their reference. 
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Figure 6. Output characteristics of the transmitter using 16-QAM OFDM signals (a) PA’s output PSD and 
(b) AM-to-AM and AM-to-PM characteristics for class AB power amplifier with 1.07 dB gain imbalance and 
5° phase imbalance. 


Table 2. Performance of proposed model in comparison with the other structures in literature. 


Ref [19] [10] [16] [17] [27] [28] This work 
Input Signal LTE OFDM 16 WCDMA WCDMA LTE 16 QAM OFDM 16 
QAM QAM 
PAPR(dB) 11.49 10 - 10.6 9.5 8 10.32 
ACPR 16 dB 21 dB 17 dB 14 dB 10 dB 11 dB 16 dB 
improvement 
Type of DPD Gegenbauer Volterra Rational Distributed Spectral Cascaded Chebyshev 
polynomial series Function polynomial Weighting DPD polynomial 
Coefficients 42 194 76 34 150 38 33 
FLOPs - - >400 333 - - 356 
Condition <30 - 117.8 68.5/48.8 - - <20 dB 
Number(dB) 
Component(s) Just PA Just PA PA+IQ PA+IQ Just PA PA+DM PA+IQ 
linearization mismatch mismatch TX mismatch 
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5. CONCLUSION 

A low complexity digital pre-distortion model was proposed to estimate and compensate for the 
defects of direct-conversion transmitters. The proposed model has compensated the nonlinearity of PA, local 
oscillator (LO) leakage, and frequency-dependent IQ imbalance of transmitter jointly, with the least number 
of DPD model coefficients and consequently, low complexity and low-cost structure for DSP implementation 
were obtained. The DPD was based on the orthogonal Chebyshev polynomial model which has been 
modified to be orthogonal for all real and complex values of the input signal and make the DPD robust to the 
input signal statistics and numerical instability. The simulation results show the promising performance of 
the proposed structure and approve its numerical stability in a fixed-point calculation environment. We 
improved the ACPR of the transmitted signal by about 16 dB for the OFDM-16 QAM excitation signals. 
Using the new orthogonal polynomial model results the conditional number parameter value reach less than 
40 decibels for the OFDM input signal. Therefore, this model can be implemented on field programmable 
gate array (FPGA) or DSP as a good candidate for the fixed-point calculation environments. 
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